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A simple one-dimensional gas-piston kinetic model gives the interaction potential between two 
colliding heavy ions. In the frame of the classical, thermodynamical approach, the colliding heavy 
ions are not submitted to friction, but produces an irreversible phenomena with cause at the differ¬ 
ence of pressure p “felt” by the nucleon gas when ions collide with the target when compared with 
the pressure that nuclear matter exert on their boundaries when in thermodynamical equilibrium, 
and offers a straightforward way to calculate interacting potentials. 
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I. INTRODUCTION 


The numerical calculations implicated in the study of 
nuclear reactions between heavy ions have an extreme 
complexity. Some times it is advisable to build simple 
models that allow to obtain concrete images of the phe¬ 
nomena based on classical conceptions and made useful 
predictions. In this perspective, D. H. E. Gross [lj pro¬ 
posed the paradigmatic piston-gas model. However, our 
interpretation differs, not sustaining a friction process. 

We propose in this Letter a classical, thermodynam¬ 
ical approach to determine the interaction potential to 
which two colliding heavy ions (with mass number A 
above 40) are submitted, and its dependency relatively 
to their nuclear temperature T. The high number of nu¬ 
cleons present inside each nucleus allows a macroscopic 
treatment. Also, the high value of the relative momen¬ 
tum ( Hk ) with k « 10 — 40 fm -1 ensures that the reduced 
wavelength tobeA = l/fc<Cl fm, which implies that 
the wave packet associated to each ion is much smaller 
that the ion radius (« 1.2A 1 / 3 fm). 

The collision reactions are irreversible and can be 
treated as a dissipative phenomena due to “friction” that 
are supposed to be occurring at the moment of impact. 
However, we introduce here a different perspective, based 
on previous work done by R. de Abreu et al. 0,13, con¬ 
sisting in the following assumption: there is no friction, 
but an irreversible phenomena with cause at the differ¬ 
ence of pressure p' “felt” by the uucleon gas when ions 
collide with a target when compared with the pressure 
that nuclear matter exert in their boundaries when in 


thermodynamical equilibrium p: 

p ' =! ’( 1 “?)’ (1) 

with r and u denoting, resp., the speeds of the piston 
and of the gas particles, supposed to collide in chaotical 
manner with the piston and in (relatively) great number. 

It is known that the variation of energy due to the 
piston is given by 

dE = -dWdiss = p'dV = p ^1 —— ^ Adr. (2) 

Notice that we are treating here the repulsive region of 
the potential when the ions are compressing each other, 
and distant from the repulsive core, r > a, otherwise we 
must consider the repulsive Coulomb potential as well. 
As a first approximation, we assume p independent from 
r, and we can immediately integrate Eq. [2] from r = a, 
with r denoting a minimal approaching distance and r a 
given distance from the center of the nucleus: 

2 t * \ dv 

1-) —dtB(T, p, m N ), (3) 

u J dt 

where A is the area of contact between the two nucleus 
and B(T,p, tun) is the pressure exerted by the nuclear 
matter, dependent on the gas temperature T, chemical 
potential p, and the nucleons mass mjv- We will assume 
that the gas of nucleons is a Fermi gas without interaction 
between its particles. The chemical potential is given by 
(see, e.g., Ref. (3|) 

o _ 2 gVrnJj 2 /'°° y/ede 

3 \Z2TT 2 h 3 Jo eT + 1 

It can be shown that 

dn _ VTgy/^m^ 2 

df ~ m 


E = ip(r) = A ^ 
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and the entropy is (S = —(dfl/dT)v,^i) 


S = VT 


g^/2jlm 
6 H 3 


3/2 
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The pressure is straightforwardly obtained: 

T 2 gV2ilm 3 ^ 2 


P = P{T = 0) + 


12H 3 


and it can be rewritten using Eq. [6] 

T S 


p=p{T = 0) + 


2 V 
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We assume that p(T = 0) ~ 0. Eq. [5] shows that the 
resulting compression due to the heavy ion collision is 
equivalent to an amount of work TS/ 2, and generation 
of entropy due to the new arrangement of particles. From 
the above, we obtain the interacting potential 


ip(r) = —A 


gT 2 


m 3 


1 - 


2r\ dr 
dt 


dt. 


( 9 ) 


Keeping with the analogy between the system piston-gas 
and the collision between heavy ions, we assume that the 
the energy of the incident ion is given by 


1 , 

£ = - Mr 2 + ip(r). 


( 10 ) 


It results 


<p(r) = - AB(T , n, m N ) j ^1 - ^ \/e- <p(r')^ dr', 


with 


tdcji n _ gT 2 V^p^n 2 

B(T, p, m N ) = 


12h 3 


( 11 ) 

( 12 ) 


We can go a step further to obtain 

ip(r ) = —AB(r — a) + — yj^AB J y/e — tp{r')dr' (13) 

where M denotes the mass of the heavy ion and a is the 
distance of minimal approximation, corresponding to the 
moment when the ions collide inelastically, dissipating 
energy and ejecting nucleons in the process. 

The mathematical resolution of Eq.[13]is complex since 
we don’t have a previous idea of what is the form of the 
interaction potential tp(r) (in fact, it is what we want to 
know). However, we may assume that the total energy of 
the heavy ion is much higher that the interaction poten¬ 
tial in the range considered (until the minimal distance), 
e tp(r) (relativistic heavy-ions collisions), in which case 
we obtain a more simplified integral equation: 

ip(r) = f(r ) + A f ip(r)dr. (14) 

J a 



FIG. 1: Ion-ion potential for reaction (a) ls O + ls O and (b) 
90 Zr + 208 Pb. The touching point distance of two ions is 
a = 6, 9 fm, </>(a) = —5.77 MeV, C = 0.26. The square dots 
are data representing the calculations obtained in the Hartree- 
Fock-Bogoliubov approximation with SkM* parameter set of 
the Skyrme force [jjj]. 


Here, 


f( r ) = -[AB(T,fj,,m N ) + y/e\(r - a), (15) 

and 

x = -~\ TT AB ( T ’ /b"bv)- (I 6 ) 

u V Me 

In the relativistic limit, the mathematical solution of 
Eq. [14] is given by 


V(r) = ^ + v(a)e-^ r ~ a l (17) 

C is a constant of integration that appears after convert¬ 
ing Eq. [13] to a first-order differential equation. In Fig. [T] 
it is shown the calculations using Eq. [l7] for the ion-ion 
potential for reaction 16 O + 16 O. The touching point 
distance of two ions is a = 6.9 fm for the set of three 
reactions shown in Table U despite the different number 
of nucleons in each nucleus, cj)(a) = —7.0 MeV, C = 0.26 
Mev/fm, and | A |= 1 fm -1 . Notice that r c = 1/ | A |= 1 
fm is the critical distance in remarkable agreement with 
the expected value @. In the Woods-Saxon spherical po¬ 
tential describing single-particle levels, r c represents the 
surface thickness (e.g., Ref. 0 )- 

The comparison is made with several calculations ob¬ 
tained in the Hartree-Fock-Bogoliuobov approximation 
with SkM* parameter set of the Skyrme force [Ej and 
our theoretical approximation of Eq. [17] with parameters 
shown in Table |T] 

The advantage of the method proposed in this Letter 
stays in the possibility to obtain an analytical form for 
the interaction potential, more easy, fast and precise to 
use. The simple physical model also points out that what 
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FIG. 2: Ion-ion potential for 90 Zr + 90 Zr. The square dots 
are data representing the calculations obtained in the Hartree- 
Fock-Bogoliubov approximation with SkM* parameter set of 
the Skyrme force [a]. 


TABLE I: Parameters for the interaction potential in nuclear 
heavy ions collisions. Potentials best fitted for A = 1 fm _1 . 


Nuclear reaction 

a (fm) 

(j>{a) (MeV) C (MeV/fm) 

16 o + 16 o 

6.9 

-7.0 

0.26 

90 Zr + 90 Zr 

6.9 

-1330.0 

0.26 

90 Zr + 208 Pb 

6.9 

-40.0 

0.26 


occurs between the colliding heavy ions is not a friction 
process, but simply an irreversible process due to the 
difference of pressures felt by the nucleon gas before and 
during the collisional process. 


II. CONCLUSION 


It was obtained the interaction potential between two 
heavy ions in the framework of a simple, classical, one¬ 
dimensional piston-gas model, where there is no friction, 
but an irreversible phenomena with cause at the differ¬ 
ence of pressure p' “felt” by the nucleon gas when ions 
collide with a target when compared with the pressure 
that nuclear matter exert in their boundaries when in 
thermodynamical equilibrium p. Our parametrized po¬ 
tentials are in good agreement with other theoretical po¬ 
tentials, particularly with the Hartree-Fock-Bogoliuobov 
approximation with SkM* parameter set of the Skyrme 
force, and with experimental data. 
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